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1. Introduction 



Scalar condensates often play an important role in cosmology. For example, the 
superluminal expansion of the early inflationary universe is driven by primordial dark 
energy, which in scalar field inflation is the energy of the zero mode inflaton condensate 
PQ. At the end of inflation the condensate <fr oscillates and then decays roughly speaking 
when the decay rate > H, where H is the Hubble rate. The decay products are 
assumed to thermalize almost instantaneously by scattering against each other, thus 
reheating the cold universe. While the inflaton condensate is decaying, the expansion 
rate of the universe depends on the varying relative magnitudes of the plasma and the 
condensate energies, which thus affect the decay and reheating processes. It turns out 
[21 EH E] that there is a maximum temperature Tmax which the plasma can attain; this 
is reached when the plasma energy is still subdominant. The reheat temperature Trh 
marks the time after which the universe is dominated by the relativistic gas of the decay 
products. 

The process of thermalization of a scalar condensate in the heat bath of its decay 
products in an expanding universe is a generic problem that is relevant not only for 
inflation but, for example, also for the decay of the squark/slepton condensate of the 
minimally supersymmetric standard model (MSSM) j^J. In thermalized plasma the 
dispersion relations of the decay products are modified by what is commonly known 
as plasma masses |.6^ (the role of plasma masses in early universe was first noted in 
[7j). These can affect the decay rate, as was recently discussed in jH] in the context 
of inflation. There the plasma masses were put in by hand into the decay rate 
to provide a kinematic suppression factor. The authors found that the inclusion of 
the plasma masses gives rise to an era of constant temperature during reheating; this 
happens because the masses of the decay products of the inflaton increase with the 
temperature and at some point the energy of the decay products would exceed the 
energy of the inflaton zero mode, rendering the decay kinematically forbidden. 

However, a thermal background induces modifications in the scalar condensate 
decay rate which are more subtle than a simple phase space suppression. If one assumes 
that the scalar condensate decays mainly into fermions, as we will do in this paper, the 
plasma has two kinds of elementary fermionic excitations, called particles and holes. 
Particles are modes which in the limit T — > yield ordinary elementary excitations, 
whereas holes are collective modes present in thermal bath only; they are called "holes" 
since these excitations correspond to the removal of an antiparticle, creating a state 
with all the same quantum numbers as a particle state, but with different energy j^j. 
Thus for instance an inflaton condensate can decay into a pair of particles or a pair of 
holes, or it can transform into a particle by absorbing a hole from the thermal bath. All 
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these processes contribute to the total thermalization rate of the inflaton condensate 
and need to be evaluated. 

In this paper we calculate the thermalization rate of a scalar condensate decaying 
into fermionic particles at finite temperature. The results could be applied to any 
situation where there is an initially dominant scalar field decaying in an expanding 
universe. Here we focus on the reheating phase of the early universe. 

2. Finite temperature thermalization rates 

Let us consider the interaction Lagrangian 



where is a some scalar field that forms a zero mode condensate in the early universe, 
and if> is a fermion. Here gy is the Yukawa coupling between the scalar and the 
fermion. In what follows we will neglect the self-couplings of as well as its couplings 
to other scalars and assume that the fermions if) are brought to thermal equilibrium 
instantaneously after they are produced by virtue of if>tf> scatterings. We also assume 
that the rest masses of the fermions are much smaller than the ambient temperature 
T, so we can take the fermionic dispersion relation to be dominated by the first order 
temperature correction m = CgT, where g is a coupling constant appropriate for ipif) 
scattering and C is some number, which we take to be 1 for simplicity. 

The thermalization rate r^(T) of a scalar condensate has been computed in the 
past (see e.g. ^U]) and is given by the imaginary part of the zero mode propagator E; 
at one- loop level this is the cut of the relevant one-loop diagrams at finite T. In general 
one can write 



where u is the energy of the given mode. £ receives contributions from all the fermionic 
modes present at finite temperature. 

When considering the thermalization of a scalar field, there are two kinds of 
processes. First, the scalar zero mode can decay directly into a pair of particles or 
a pair of holes. The energies of the produced particles are then constrained by 



where p denotes particle and h hole. Mostly we will be interested in a cold condensate 
dominating the universe, such as the inflaton or the MSSM flat direction condensate. 
Hence we may ignore thermal corrections to the scalar mass and take u = M^. 

Second, if the mass of the scalar is small enough (or temperature large enough), the 
possibility opens for a hole in the heat bath to absorb the scalar to produce a particle. 



£ = gycpipip , 



(1) 




(2) 



= 2ui(k) ; i =p,h , 



(3) 
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This could be called thermal fragmentation of the condensate. In this case the energies 
are constrained by 

M <t> + u h (k)=u) p (k). (4) 

As will be discussed below, it turns out that for thermal fragmentation to occur, 
the temperature of the thermal bath has to become larger than the scalar rest mass. 
However, since the scalar-fermion (heat bath) coupling, gy, is typically much smaller 
than the coupling g relevant for ipip scattering, we can ignore the thermal mass of the 
scalar. 

The dispersion relations for fermionic positive energy hole and particle excitations 
with thermal masses m = gT can be written as [HI HH] 
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where u 
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uj/T, k = k/T. The dispersion relation for holes can be also expressed as jH] 

~P k 
— + - — 

9 uj h + , 

this form turns out to be more convenient in numerical calculations. The relations Eqs. 
©-©) together with the energy conditions for absorption and decay, can be solved 
numerically to obtain u>h, uj p and k as functions of temperature, and the results can be 
used to find T^(T). The imaginary parts needed for the decay rates were computed in 
HU and read (t = T/Mj) 
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for the case of absorption, whereas 
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(9) 



for the case of decay; here gy is the Yukawa coupling between the scalar and the fermion 
and 

1 



(10) 



exp^p) + 1 
is the distribution function. 

The total scalar condensate thermalization rate Ta(T,u>) = r abs + T dec is a 
complicated function depending on the coupling constant g. However, varying g and 
fitting the numerical results with simple functional forms we find that the thermalization 
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rates are well described by 

T dec (t,9) = ^(4 - 6gt) (ae- b ^ 2 + c) (l - ^^j) , (H) 

r o6s (t, g) = ({D 2 g 2 + Dig + D )t + ^ , (12) 

where a = 0.34; b = 5.6546; c = -0.09 and D 2 = 3.746; D 1 = -0.7058; £> = 
11.1278; E = 1.8015. 

Since the direct decay channel is the only possible one when T is small enough, 
while the absorption channel opens only when T is sufficiently large, there is a range of 
temperatures where the both reactions are kinematically forbidden if the scalar couples 
only to a single fermion. In reality the condensate is more than likely to couple to many 
species of fermions. For instance, the MSSM flat direction condensate typically couples 
to meny quark and/or lepton generations. Likewise, the inflaton could well couple 
to several fermionic species. If a condensate couples to more than one species with 
different thermal masses, it is possible to obtain a piecewise continuous non-zero decay 
rate. This is because the end point of the direct decay channel and the starting point 
of the absorption channel depend on the coupling constant: if the scalar condensate 
decays directly to a fermion with a thermal mass proportional to a coupling constant 
gi and absorbs a fermion with a thermal mass proportional to a coupling constant g 2 , 
then in order for the temperature gap in the thermalization rate to vanish the ratio of 
couplings must satisfy the relation 

— < 0.25 , (13) 
9i 

which is easily realized in actual models. The Yukawa couplings may of course also vary, 
but here we keep all the Yukawas equal for simplicity. 

To demonstrate the interplay of the decay and absorption channels in case of more 
than one fermion species, we have calculated the decay rates for a scalar coupled to two 
fermions with thermal masses proportional to coupling constants g% = 1, g 2 = 0.1 and 
9i — 1) 92 — 0.2; the results are plotted in Fig. [TJ 



3. Solving the Boltzmann equations 

The energy partitioning between the scalar condensate and the relativistic heat bath and 
its time evolution in the expanding universe are described by the Boltzmann equations, 

p + 3Hp (f> + r^,(T)p0 = , 

p R + AH PR - T^T) P4> = , (14) 
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Figure 1. The complete thermalization rates for a scalar zero mode decaying into 
two fermionic particles with plasma masses proportional to g\ — 1 and <?2 = 0.1 (left), 
and g± = 1 and 172 = 0.2 (right). The low temperature part of the plot corresponds 
to the direct decay channels, while the zigzag patterns are due to absorption channels 
opening first to particles with the largest thermal mass, then to particles with a smaller 
thermal mass. 



where and pR are respectively the inflaton and the radiation energy densities. Using 
the dimensionless quantities 

$ = P4> M^a 3 ; R = p R a 4 (15) 

as was done in e.g. (HJ, and defining a new variable x = aM^, Eq. (j!4j) can be written 
in a form more suitable for numerical calculations (' = d/dx), 

$' = _ ,/X^r,m _ x cf> (16) 
r> = #^r,(T) * 2 9. (17) 

Assuming that at the beginning the condensate energy density dominates, the initial 
conditions at x = xj are given by R(xi) = and $(x/) = $/, where $7 is obtained from 
the Friedmann equation, 



8tt Ml MT 



The temperature of relativistic gas in an expanding universe is given by 
T(x) ( 30 \ 1/4 R 1 ^ 



(19) 



where g* ~ 100 counts the effective degrees of freedom. The time evolution of 
temperature as a function of the scale parameter x, computed with T^T) of Eq. (j2J) 
is plotted in Fig. E] together with the results presented in [B]. We have calculated the 
temperatures for two sets of initial values following the cases II and III discussed in 
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Figure 2. Solutions for Boltzmann equations with two different sets of initial values: 
g\ = 6 X 1CT 13 , M = 10 9 GeV, V 1 ^ = 10 14 GeV (left) and g\ = 3 x 10~ 9 , 
M = 2 x 10 7 GeV, V" 1 / 4 = 8 x 10 n GeV (right). The curves are: results obtained 
by Kolb et al. in [5j (thin solid); the case where thermal corrections are neglected, i.e. 

= 9y^4> (thin dashed); the case where thermal corrections are included with gx = 1 
and 92 = 0.1 (thick dashed); and the case where thermal corrections are included with 
gi = 1 and gi = 0.2 (thick solid). 

first with the scalar mass falling between the reheat temperature and the maximum 
temperature reached by the universe during the reheating, Trh < < Tmax, with 
the reasonable numerical values g\ = 6 x 1CT 13 , = 10 9 GeV, V 1 ^ 4 = 10 14 GeV, 
where V is the value of scalar potential at the end of inflation; when the scalar mass is 
smaller than the reheat temperature, < T RHl we adopt the values g\ = 3 x 10" 9 , 

= 2 x 10 7 GeV, and V 1 ^ = 8 x 10 n GeV. 

As can be seen in Fig. El an era of a constant temperature which was found by 
a partial consideration of the thermal effects relevant for the decay of the inflaton 
condensate jH], actually vanishes when one takes into account all the interactions 
between the scalar and the thermal bath at the one-loop level. The inclusion of the 
absorption decay channel and the possibility to decay into more than one type of fermion 
allow for the reactions to continue at temperatures greater than the inflaton rest mass. 
This effect is due to the existence of several fermionic species coupling to the inflaton; 
in the case of one species only, the result would be almost identical to the case of Kolb 
et al. jH], with only minor changes due to the hole excitations. As a consequence 
the thermal history of the heat bath turns out to be qualitatively similar to the case 
when the thermal masses are ignored. At the end of inflation the rapid decays of 
inflaton raise the temperature to some maximum temperature, after which most of the 
time the temperature scales as T oc a" 3 / 8 ; however, there are short periods during 
which the temperature drops faster because of absorption channels reaching the level 
of fastest direct decay rates before closing with decreasing temperature. If the thermal 
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bath contained also bosons that couple to the inflaton, the thermalization rate would 
be modified accordingly. The details would again depend on the coupling strengths. 
Clearly, the thermalization rate and the temperature of the resulting heat bath very 
much depend not only on the details of the inflaton model itself but also on the reaction 
rates of the decay products, which determine their plasma masses. Such considerations 
might also be important for e.g. thermal leptogenesis [TTj . 
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